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We study the effect of gravitational focusing of the earth on dark matter. We find that the effect 
can produce a detectable diurnal modulation in the dark matter signal for part of the parameter 
space which for high dark matter masses is larger than the diurnal modulation induced by the 
fluctuations in the flux of dark matter particles due to the rotation of the earth around its own axis. 

The two sources of diurnal modulation have different phases and can be distinguished from each 
other. We demonstrate that the diurnal modulation can potentially check the self-consistency of 
experiments that observe annual modulated signals that can be attributed to dark matter. Failing 
to discover a daily varying signal can result conclusively to the falsification of the hypothesis that 
the annual modulation is due to dark matter. We also suggest that null result experiments should 
check for a daily modulation of their rejected background signal with specific phases. A potential 
discovery could mean that dark matter collisions have been vetoed out. 


I. INTRODUCTION 

Direct detection of dark matter (DM) particles requires 
sophisticated ways of distinguishing a potential DM sig¬ 
nal from noise that is due to other not so exotic particles 
such as neutrons, muons, photons etc. Most of direct 
detection experiments rely on vetoing signals that are 
consistent with the background. Therefore the correct 
and precise rejection of the background signal can make 
the difference between discovery of DM or not. Needless 
to say that although experiments have made enormous 
progress in understanding the respective backgrounds, 
one cannot be absolutely sure that the appropriate sig¬ 
nal has been vetoed. However, there are other ways to 
circumvent the background problem. Experiments like 
DAMA rely on a different way of probing DM. Instead of 
rejecting background signals, experiments like these are 
looking for an annually modulated signal which could be 
consistent with DM. The earth is moving with respect to 
the reference frame of the DM halo of our galaxy with 
a velocity that is mainly the velocity of the solar system 
superimposed by a much smaller component, i.e. the 
velocity of the earth around the sun. A variation in the 
velocity of the earth with respect to the DM halo leads to 
a variation of the flux of DM particles that reach an earth 
based detector, thus causing the aforementioned annual 
modulation of the signal. If no background component 
has a similar annual modulation, the annual modulation 
is solely attributed to DM. In principle it does not matter 
if someone is agnostic regarding the experimental back¬ 
ground, since the existence of an annually modulated sig¬ 
nal with a phase expected from the simplest models of 
the dark halo would be a strong hint that this is only a 
DM effect. Based on this principle, DAMA has observed 
an annual modulated signal that is consistent with the 
existence of DM [1]. 
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However there can be more sources of modulation in 
the DM signal. An obvious one is the angular velocity 
of the earth around its own axis. The earth is rotating 
around its own axis and therefore, the net velocity of the 
earth with respect to the DM halo fluctuates daily [2]. 
Given that the earth is moving with an average speed of 
vq = 220km/s (in the local standard of rest), the fluctua¬ 
tion due to the angular velocity of the earth is tiny, since 
the speed of a point on the equator is ^ .5km/s. If the 
detected signal from DAMA is due to DM-detector col¬ 
lisions, one should at some point be able to observe also 
the diurnal modulation of the DM signal. This is an im¬ 
portant self-consistency check because if DAMA observes 
an annual modulation and does not observe at some point 
a diurnal one, this is inconsistent and an alternative ex¬ 
planation of the annual modulated signal must be found. 
This is true for all experiments that base their detection 
technique on annual modulated signals. In fact recently, 
the DAMA collaboration studied the aforementioned ef¬ 
fect concluding that it has not reached yet the required 
exposure to observe the diurnal modulation of the DM 
signal due to the angular velocity of the earth [3]. An¬ 
other source of diurnal modulation relevant for low DM 
masses is due to stopping via DM interactions with un¬ 
derground nuclei [4] . Depending on the distance that DM 
particles travel underground in order to reach the detec¬ 
tor, they lose different amounts of energy due to interac¬ 
tions with atoms and nuclei. Due to the rotation of the 
earth, this creates a small amount of diurnal modulation 
signal. This can be particularly important in the case 
where DM is in the form of Strongly Interacting Massive 
Particles [5, 6] or mirror dark matter [7, 8]. Moreover, di¬ 
urnal modulated DM signals can emerge from channelling 
effects in the detectors, since DM-nuclei scatterings from 
specific directions might have no quenching [9, 10]. Addi¬ 
tionally, there can be a modulation in the DM signal due 
to gravitational focusing [11]. Massive objects can work 
as a gravitational lens, changing the DM phase space 
density around them. This effect has been studied in the 
context of the sun acting as a lens, thus producing an 
annual modulated signal since the position of the earth 
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with respect to the DM flux and the sun changes during 
the year [12, 13]. One should note here that although 
this amount of annual modulation is much smaller com¬ 
pared to the one due to the velocity of the earth revolving 
around the sun, it has a different phase and therefore it 
could potentially be distinguishable. Similarly, monthly 
modulated signals are also probable due to the focusing 
effect of the moon as well as the monthly varying velocity 
of the earth around the center of mass of the earth-moon 
system [14]. 



t = 18 hr 


FIG. 1: This cartoon illustrates the two diurnal modulation 
signals of a detector rotating in the dark matter wind (dashed 
lines). We choose t = 0 hr to be the time where the detector 
is maximally aligned against the dark matter wind. The red 
points are the extrema of the angular velocity modulation. 
The blue points at t = 0 hr and t = 12 hr are the minimum 
and maximum of the focusing modulation. 

In this paper we study the effect of gravitational fo¬ 
cusing produced by the earth itself acting as a lens. This 
produces a diurnal modulated signal because as the earth 
rotates around itself, there is a different amount of fo¬ 
cusing depending on the relative position of the detector 
with respect to the center of the earth and the direction of 
the DM wind (i.e. the direction where the velocity of the 
earth with respect to the reference system of the DM halo 
is maximum, see Fig. 1). We find that the effect can be 
significant when the DM mass is large. Assuming spin- 
independent contact interactions and elastic scattering 
the DAMA annual modulated signal has roughly two re¬ 
gions in the cross section-DM mass parameter space that 
are consistent with DM, i.e. a low mass region around 
10 GeV and a heavier one ^ 50 GeV. We find that the 
diurnal modulation due to earth gravitational focusing is 
larger than that produced by the angular velocity of the 


earth for part of the DM mass parameter space and there¬ 
fore it should not be ignored. As we will demonstrate 
the diurnal modulated signal due to focusing is currently 
“scratching” the allowed parameter space of DAMA in 
the large DM mass region. We should also emphasize 
that the two sources of diurnal modulation produce sig¬ 
nals with different phases, therefore they can easily be 
identified once sufficient exposure has been reached. 

In the next section we review the effect of gravitational 
focusing providing the relevant formulas. In section III, 
we calculate the diurnal modulation for gravitational fo¬ 
cusing and due to the angular velocity of the earth, pro¬ 
viding the amount of expected counts. We conclude in 
section IV. 


II. GRAVITATIONAL FOCUSING 

The distortion of a distribution of particles due to the 
presence of a massive object has been studied analyti¬ 
cally [15-19]. Here we are going to follow the derivation 
of [18] the validity of which has been verified by numeri¬ 
cal simulations [19]. 

Let us assume that DM follows a velocity distribution 
/(zToo) in the absence of any focusing effect. The distribu¬ 
tion is obviously independent of the position. Liouville’s 
theorem states that the phase space density is conserved 
along particle trajectories. This means that the distribu¬ 
tion of DM in the presence of a massive object (e.g the 
earth) should be 


fe{r,v) = f{vao{r,v)), (1) 


where Vao (u v) is the velocity far from the gravitational 
source of the trajectory that passes by the position r with 
velocity v. We briefly review how one obtains Vca follow¬ 
ing [18]. In order to obtain this velocity one should use 
the known constants of motion along the trajectory of the 
particle from a far away position where the gravitational 
field of the earth is negligible to the point where it enters 
the earth. The constants of motion are the energy, the 
angular momentum and the Laplace-Runge-Lenz vector 
per unit DM particle mass: 
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where is the mass of the earth. Setting the values of 
E, I and A at asymptotically far away distance equal to 
the corresponding values at distance f from the center of 
the earth, one finds [18] 
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where P = — (r • z7)^, a = GMe/v"^, Voo is given by Using Eq. (1) we get 
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We use a truncated Maxwell-Boltzmann DM distribution 
with Voo < Vesc where v^sc is the escape velocity from the 
galaxy. At asymptotically far away distance from the 
earth it takes the form 


feir,v) = 
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where A^esc = erf (uesc/'^o) - 2uesc exp(-u2s(./ug)/(i/7ruo) 
and Ve is the velocity of the earth with respect to the rest 
frame of the DM halo. The velocity vq is given by 
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where S is the angle between and the DM velocity v 
and /3 is the angle between Ve and the position r of the 
observer relative to the center of the earth. We will take 
r to be a point on the earth’s surface. After a particle 
has reached the surface of the earth, it will move under¬ 
ground until it reaches the detector. However for this 
part of the trip, the above equations do not hold. This is 
because the derivation of distribution (6) describes grav¬ 
itational focusing around a point mass Mg. After the 
particle has entered the earth, one cannot treat the earth 
as a point mass. Additionally, the Laplace-Runge-Lenz 
vector of Eq. (2) is no longer conserved because inside 
the earth the gravitational force does not scale anymore 
as ^ r~^ which is a necessary condition for Eq. (2) to 
be conserved. For the remaining part of the DM parti¬ 
cle trajectory inside the earth we assume a constant mass 
density which corresponds to a gravitational force scaling 
as ^ r. For such a force there is a new Laplace-Runge- 
Lenz vector that is conserved. In the appendix, by apply¬ 
ing conservation of energy, angular momentum and the 
new Laplace-Runge-Lenz vector, we demonstrate that in 
the undeground part of the particle’s trip, no significant 
enhancement in the focusing takes place. Therefore, it is 
an excellent approximation to track the trajectory of the 
particle from asymptotically far distances to the surface 
of the earth and then assume that the particle travels in 
a straight line until it reaches the detector. In any case 
our approximation lies on the conservative side i.e. the 
actual focusing can only be larger than that calculated 


within our approximation since we ignore the focusing 
effect for the part of the particle’s trajectory inside the 
earth. As we argue in the appendix, this makes a negli¬ 
gible difference in practice. 

In order to compute how gravitational focusing of the 
earth can induce a diurnal modulation of the DM signal 
we follow the convention of [4] and introduce 

h = xcos9i cosujt + {jcosOl sinwt ± zsin^;, (8) 

De = xsina-I-zcosa, (9) 

where n is the unit vector with direction from the center 
of the earth to the detector, Ve is the velocity of the earth 
in the rest frame of the galaxy, 9i is the latitude of the 
detector, and a the angle between Vf. and the angular 
velocity w of the earth. We have chosen a galactic rest 
frame where the z-axis is along the north-south pole, and 
Ve lies along the x — z plane. It is understood that the 
± corresponds to the north and south hemisphere. We 
have chosen t = 0 the time where Ve and h align as much 
as possible i.e. fi is along the x — z plane. Recall that for 
the evaluation of Eq. (7) we need to know both angles 
6 and (3. By writing the velocity of the DM particle in 
spherical coordinates 

V = ^(aisin^cos/!) -I- y sin 0 sin </> -|- zcos0), (10) 

we can express the angle 5 as 

cos 5 = V ■ Ve = sin a sin 0 cos 0 -I- cos a cos 0. (11) 
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In order to estimate the angle /3 (defined as cos /3 = r-Ve), 
we need to associate the point where the particle crosses 
the surface of the earth with coordinates f with other 
known angles. From the geometry of the problem and 
upon assuming that particles move on straight lines after 
they have crossed the surface of the earth, it is easy to 
see that 


r = 



Rffi — (-D ^ 



( 12 ) 


where i ?0 is the radius of the earth, Id the depth of 
the detector, and t the length of the trajectory (i.e. the 
length of the straight line) from the point where the par¬ 
ticle enters the earth to the detector. This length is given 
by 


where is the DM mass, = 0.3 GeV/cm^ is the DM 
density, Vesc = 550 km/s is the galactic escape velocity 
and is the number of nuclei targets of type i in the 
detector. For example in the case of DAMA, the detector 
is made of Nal where iV 7 ’,Na = Nt.i = A/tot/(wNa + wi) 
(Mtot being the total mass of the detector). /e(r,F) is 
given in Eq. ( 6 ), and da/dEji 

is the differential DM-nucleus cross section, rriN and p, 
being the nucleus mass and the reduced mass of DM and 
nucleus respectively. 

The diurnal modulation caused by the angular velocity 
of the earth spinning around its own axis enters through 
the time dependence of Ve in Eq. (5). Ve is 

Ve{t) =Vq+ Vrevit) + F-ot (f) , ( 18 ) 


£ = {R(^ — £d) costp (13) 

+ ^(i?® - £d? cos 2 F - {£l - 2RM, 

where t/j is the angle between the velocity of the particle 
and n. From Eqs. ( 8 ) and (10) we get 

cos Ip = V ■ n = cos 9i cos uit sin 9 cos (p 
+ cos 9i sin wt sin 6* sin ^ ± sinF cos0. (14) 

Using Eq. (12) we get 

P - - ^ ^ , R® - (-D ^ ^ 

cos p = r ■ Ve = cos 0 H- - - n ■ Ve- (15) 

it® it® 

From Eqs. ( 8 ) and (9) is trivial to deduce n • Ve and the 
final expression for cos j3 becomes 

cos /3 = 

cos <5 + ^(sina cos 0 cos wt ± cosasinF). 

it® it® 

(16) 


where F® is the velocity of the sun in the galactic frame, 
Frev(F i® earth’s velocity as it revolves annually around 
the sun, and Frot (i) is the rotational velocity of the earth 
around its own axis. Since Frev(l) varies slowly and is 
small compared to F® the term relevant on daily time 
scales is Frot(i)- We will take Vrev{t) — 0 since it will only 
contribute to the constant part of the velocity during any 
particular daily cycle. In absolute values we write: 


Ve{t) C±Vq + Urot COS [w(t - Uot)] , (19) 

where Vq = 232±50 km/s, w = 27r/24 hr“^ and tj-ot = 18 
hr. Note that our t = 0 corresponds to the time where the 
detector aligns maximally against the DM wind. VJ-ot = 
wi?® cos 9i cosK = 0.23 km/s which is the speed due to 
the rotation of the earth at latitude 9i = 42 degrees, k 
being the angle between F^ot and F® [3]. 

We assume that scattering on detector nuclei is de¬ 
scribed by a contact interaction 
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One can see from Eq. (7) that vc depends on /? which 
depends on time. This is what creates the diurnal mod¬ 
ulation in the DM signal. 

III. DIURNAL MODULATION 


where we provide the cross section in terms of the normal¬ 
ized DM-nucleon cross section tjp. pp is the DM-nucleon 
reduced mass, A is the number of nucleons in the nucleus 
and F(Eji) is nuclear form factor. If one neglects gravi¬ 
tational focusing for the moment, the velocity integral of 
Eq. (17) can be calculated analytically 


The differential rate of dark matter events in a detector 
is given by 
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The result of this integral is [21]: 
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FIG. 2: The solid (dashed) lines correspond to a modulation amplitude of 10“® (10“^) cpd/kg/keVee in the 2-6 keV bin 
for DAMA. The red lines correspond to the modulation induced by the angular velocity of the earth, whereas the blue ones 
correspond to gravitational focusing. The green contours are the region favoured by DAMA at 3cr and 5a. 


where x = Vmin/vo, y = Ve{t)/vQ and z = {v^sc + Ve)/vo. 

Although our formalism is generic and can be valid for 
any experiment, we are going to focus here on DAMA 
since it is the experiment with a confirmed annual mod¬ 
ulation signal with more than 9a confidence when con¬ 
sidering the full exposure of 1.33 ton x yr. The DAMA 
collaboration has looked for diurnal modulated signals 
in [3] where an upper bound for the amplitude of diur¬ 
nal modulation has been imposed, albeit only consider¬ 
ing the effect of the angular velocity of the earth and 
not the gravitational focusing. In this paper DAMA has 
presented the experimental diurnal residual rate of the 
single-hit scintillation events, in different energy bins as 
a function of the hour of either a solar or a sidereal day. 
The cumulative exposure analysed in [3] is 1.04 ton x yr. 
In order to compare with the registered counts of DAMA 
we use 

where G{E'^,Er) = eyL^\-{E'j^- qE^Y/2al]/V^ae is a 
convolution function that takes into account the energy 
resolution of the experiment, q is the quenching factor 
of the scattering nucleus (i.e. 0.3 for Na and 0.09 for 
I) and (Te characterizes the energy resolution of the ex¬ 
periment [20]. The total rate in the fc-th energy bin is 
then 

R, = / ;^dE'R, (24) 

JEk dEj^ 

and the modulation amplitude in cpd/kg/keVee is 

^ _ 1 maxt - mint Efc 


Fig. 1 illustrates the fact that the maximal event rates 
for gravitational focusing and the rotation of the earth 
occur at different times. For gravitational focusing max¬ 
imum and minimum rates occur at t = 12 hr and t = 0 
respectively. For the angular velocity modulation ampli¬ 
tude in the 2-6 keVee energy bin, if DM mass is below 
152 GeV, the maximum and minimum occur at t = 18 
hr and t = 6 hr respectively simply because when the ro¬ 
tational velocity of the earth aligns maximally with Fq, 
the DM flux maximizes. However, for DM masses above 
152 GeV the maximum occurs at t = 6 hr and the mini¬ 
mum at t = 18 hr. The reason for this phase shift is that 
although the overall DM flux is maximum at < = 18 hr, 
Eq. (17) is dominated by low velocities. 

Fig. 2 shows the DM-nucleon cross section vs DM mass 
parameter space assuming the distribution of Eq. (6) and 
cross section of Eq. (20). The lines of Fig. 2 show what 
part of the DM-nucleon cross section vs DM mass leads 
to diurnal modulation rates (in the 2 to 6 keVee bin of 
DAMA) of 10“^ cpd/kg/keVee (solid lines) and of 10“^ 
cpd/kg/keVee (dashed lines) for either gravitational fo¬ 
cusing (blue lines) or angular velocity (red lines). The 
pole in the red curves corresponds to the flip in phase 
for the angular rotation effect where the modulation am¬ 
plitude vanishes and an infinite cross section is needed 
to provide a finite modulation rate. For a sense of scale 
DAMA/Libra phase-1 has already excluded modulations 
larger than 1.2 x 10“^ cpd/kg/keVee at 90% C.L. [3] (al¬ 
though this limit is only on the angular velocity effect 
with a constant phase of tmax = 18 hr). One can see 
that gravitational focusing becomes the dominant source 
of modulation at large dark matter masses. In our plot 
we have also included the DAMA favoured region based 
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FIG. 3: The predicted modulation amplitude of daily modu¬ 
lation for the angular velocity effect (red) and gravitational 
focusing (blue), calculated at dp = 10“^° cm^. The amplitude 
for an arbitrary cross section is obtained by rescaling with a 
factor Up jap. 


FIG. 5: Respective times where the maximum (x) and the 
minimum (•) number of counts occur as a function of the DM 
mass. The dashed line f = 12 hr is the time where gravita¬ 
tional focusing maximizes (if considered alone). 



FIG. 4: The predicted modulation amplitude of daily modu¬ 
lation for a Xenon detector with perfect energy resolution, i.e. 
G(F;)j,£'h) of Eq. (23) is a delta function 5{Er - E'r). We 
show the amplitude due to angular velocity (red) and gravi¬ 
tational focusing (blue), calculated at dp = 10“^'^ cm^ in the 
10-20 keV recoil energy bin. The amplitude for an arbitrary 
cross section is obtained by rescaling with a factor Up/dp. 


on the observation of the annual modulation. 

In Fig. 3 we show the predicted modulation amplitude 
in cpd/kg/keVee in the 2-6 keVee bin of a Nal detec¬ 
tor at a 42 degree latitude (the latitude of Gran Sasso), 
where a contact type interaction has been assumed. Al¬ 
though we present our results for a standard value of 
dp = 10“'^° cm^, the modulation amplitude can be ob¬ 
tained for an arbitrary cross section simply by rescaling 
the number of counts in the plot with a factor cTp/dp. At 
small DM masses (< 100 GeV), earth’s angular velocity 
is the dominant source of diurnal modulation. At large 
masses (> 100 GeV) the gravitational focusing causes a 
higher diurnal modulation rate. We find that the overall 
interference is constructive in the region where gravita¬ 
tional focusing is sizeable. In Fig. 4 we present the same 
results assuming a Xenon target. 

As we have already mentioned the two sources of di¬ 
urnal modulation we examined have generally different 
amplitudes and different phases. This is important be¬ 
cause it is possible to identify the two signals. In the cases 


where both signals are comparable, due to the difference 
in the phase, they can interfere either constructively or 
destructively. In this case the times where the maximum 
and minimum number of counts occur will not coincide 
with either one of the two individual ones (i.e. gravi¬ 
tational focusing and angular velocity). It will rather 
occur in different times according to how the two signals 
interfere. The times during the day where maximum and 
minimum number of counts are observed do not depend 
on the cross section but only on the DM mass and on 
the energy. In Fig. 5 we have plotted the time at which 
the maximal and minimal rates are observed during the 
diurnal cycle. At small DM masses the angular velocity 
effect is dominant so the maximum rate takes place at 
^max = 18 hr. At a DM mass of ~ 100 GeV gravitational 
focusing becomes the dominant effect and the maximum 
rate is found at tmax = 12 hr. Finally when the phase of 
the angular velocity effect shifts at = 152 GeV, the 
combined tmax settles at around 10 hr. The minimal rate 
occurs approximately 12 hr after the maximal. 

In addition to Fig. 3 and 4 we show in Fig. 6 the pa¬ 
rameter 


Ar?(Wmin) = ^ 


max ri{vynin, t) - mm ri{vynin, t) 


(26) 


where r] is the velocity integral in Eq. (21) and A?] 
corresponds to the difference between the maximum 
and minimum value of rj in the course of a day. For 
the modulation induced by earth’s rotation rj has the 
analytical expression of Eq. (22), while for focusing 
we numerically integrate Eq. (21) with the distribu¬ 
tion of Eq. (6). The modulation of rj is more model 
independent than the rate, since it is not integrated 
over a specific energy range and only assumes a de¬ 
pendence du/dEfi oc v~‘^ and a Maxwell-Boltzmann 
velocity distribution. For the rotation effect we show 
A77rot('Cmin) = [7?(umin, 18hr) - 77(z;inin, 6 hr)]/2 and for 
the gravitational focusing effect we show Arjfodvmin) = 
[ryCumin; 12 hr) — 7](umin, 0 hr)]/2. The cusp point in 
Fig. 2 at = 152 GeV corresponds to the Umin where 
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FIG. 6: The modulation of the velocity integral r; for the two 
effects. For low Umin the focusing effect is larger. Around 
Umin ~ 200 km/s the rotation effect changes phase and be¬ 
comes the most important effect when Umin ^ 240 km/s. 


A77rot(t^min) crosses zero in Fig. 6. This can be checked 
by taking a 6 keV recoil on iodine which corresponds to 
Umin ~ 208 km/s. 


IV. CONCLUSIONS 

In this paper we studied two sources of diurnal modu¬ 
lation in the dark matter signal. The first is due to the 
gravitational focusing of the earth, i.e. the earth acts 
as a gravitational lens that changes the local dark mat¬ 
ter phase space density. Due to the rotation of the earth 
around its axis, an earth based detector that is not placed 
at the poles of the earth will experience a daily variation 
in the amount of dark matter focusing. Additionally, the 
angular velocity of earth due to its revolution around its 
own axis causes a variation of the total velocity of the de¬ 
tector with respect to the rest frame of the dark matter 
halo, thus creating a daily fluctuation of the dark mat¬ 
ter signal. As we showed in this paper the phases of the 
two diurnal variations are distinct. For gravitational fo¬ 
cusing the minimum signal takes place when the detector 
alignes as much as possible against the dark matter wind, 
while the maximum occurs 12 hours later. On the con¬ 
trary the diurnal modulation due to the rotation velocity 
of the earth around its axis becomes maximum when the 
the rotation velocity aligns maximally against the dark 
matter wind. In the latter case for heavy enough dark 
matter masses, the phase shifts by tt as we explained in 
the text. We found that for low dark matter masses the 
“angular velocity” source of modulation is larger but the 
picture is reversed above a certain mass. 

We have provided the expected number of counts as 
well as the exact times where the maximum and the 
minimum signal should be observed for any mass and 
cross section. This could be quite important for two 
reasons: Experiments like DAMA that observe annual 
modulated signals and therefore are compatible to dark 
matter signals should pay attention to the diurnal 
signal. Failing to find it would mean that there is no 
self-consistency and the annual modulation is not caused 


by dark matter. It is an important self-consistency 
check. On the other hand, we suggest null result 
experiments like Xenon and LUX to check their rejected 
background for signs of diurnal modulation. Such a 
potential discovery would mean that dark matter signals 
might resemble for example electron events and have 
been vetoed until now. 
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APPENDIX 

We examine here how good is the approximation we have 
made in the study of gravitational focusing by assuming 
that particles follow straight lines after they have crossed 
the surface of the earth towards the underground detec¬ 
tor. In other words, we ignore the focusing effect for the 
part of the trajectory inside the earth. In order to figure 
out how good the approximation is let us assume that 
the earth has uniform density. For a particle entering 
the earth the gravitational force from the usual 1/r^ de¬ 
pendence becomes F = —{FK/Z)Gpmr where p and m 
are the uniform density of the earth and the mass of the 
dark matter particle. For this linear in distance force, 
the Laplace-Runge-Lenz vector becomes [22] 

Ae= . „ „ ^ {px£-b(TOra;R-mr£;)f}, 

yrar^uiB — mr^F + 

(27) 

where B = ui = y/k/m, k = FnCpml'i 

and E and L are the kinetic energy and angular momen¬ 
tum of the particle. Af, can be rewritten in the form 

Af. = Civ X (f X v) + C 2 r, (28) 

where Ci = m?v^r(rnr^u}B — mr^E -\- and C 2 = 

(mrujB — mrE){mr^ujB — mr^E + After using 

well known vector identities, Af. finally becomes 

Aff = {Ci+C2)f-Cx{v-f)v. (29) 

Let us consider a particle crossing the surface of the earth 
at point 1 following a trajectory that ends in the detector 
at point 2. The conservation of Ag at the two points gives 

Fifi - Ai(Di • fi)vi = r 2 h - A 2 (D • n)v, (30) 

where Fi 2 is Ci + C 2 evaluated at points 1 and 2 re¬ 
spectively. Note that we have dropped the index 2 from 
the vectors and that by definition h = f 2 - Ai _2 is Gi 
evaluated at the corresponding points. Conservation of 
angular momentum in points 1 and 2 give 

viRn^fi X vi = vicr X V, 


(31) 








(33) 


where ^c = R® ~ ^-D is the distance of the detector from 
the center of the earth. The above equation gives 

sin'll = —^sin'0, (32) 

viR(s 

where ijji {ip) is the angle between fi and Vi (n and v). 
Since we want to determine how large is the effect of fo¬ 
cusing in the trajectory during the flight of the particle 
inside the earth, let us take a trajectory where the ve¬ 
locity of the particle does not align with the direction 
of the detector fi. In this case vi and fi define a plane. 
Since the gravitational force is central, the trajectory of 
the particle will remain on this plane. Therefore we can 


express the direction of vi as 

vi = eih + t2V, 

where ei .2 are coeffcients to be determined. A first equa¬ 
tion regarding ei _2 is obtained by demanding vi ■ vi = 1. 
This leads to 

^ 1+^2 + 2eie2 cosip = 1. (34) 

We can rewrite now Eq. (30) by substituting fi by that 
of Eq.(12) (since this is entry point of particle into the 
earth) and vi by Eq. (33) keeping in mind that vi ■ fi = 
cos ipi and v ■ fi = cos ip 


^r 2 -I-eiAi cos'i/'i - Ti^^ h = ^A 2 COS'!/'- 


(35) 


If we multiply both sides for example with v we get 


r 2 -I- eiAi cos V’l - hi 


Ra 


COS Ip =[ A 2 COS Ip + ri — -£2 Ai COS Ipi 

Rm 


(36) 


Eqs. (33) and (36) form a system that can be solved in 
terms of ci and £ 2 . One can use Eq. (32) to determine 
cos'll in terms of cosip given by Eq. (14). Eor a variety 
of different angles of entry for the projectile particle, we 
have numerically solved Eqs. (33) and (36) and found 


that in all cases £2 ~ 1 and £1 ^ 0. This means that 
Di ~ D and therefore the approximation we have made 
that the particle moves on a straight trajectory for the 
part of the trip which takes place inside the earth is valid. 
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